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Abstract: We study parity violation in graviton non-gaussianity generated during in- 
flation. We develop a useful formalism to calculate graviton non-gaussianity. Using this 
formalism, we explicitly calculate the parity violating part of the bispectrum for primor- 
dial gravitational waves in the exact de Sitter spacetime and prove that no parity violation 
appears in the non-gaussianity. We also extend the analysis to slow-roll inflation and find 
that the parity violation of the bispectrum is proportional to the slow-roll parameter. We 
argue that parity violating non-gaussianity can be tested by the CMB. Our results are also 
useful for calculating three-point function of the stress tensor in the non-conformal field 
theory through the gravity/field theory correspondence. 
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1. Introduction 

The origin of the chirality of the matter sector is one of the most basic problems in parti- 
cle physics. The associated violation of parity (P) and charge conjugation (C) invariance, 
however, may not be so fundamental because such a violation can be easily realized by 
introducing a P violating hidden background field or alternatively, the invariance can be 
recovered by considering a mirror image in the hidden sector. In contrast, CP violation is 
much more profound because it implies the violation of time reversal (T) invariance under 
the assumption of TCP invariance. When this T violation is transmitted to the gravity 
sector, parity violation would occur in the gravity sector because it is legitimate to require 
the C and TCP invariance of the gravity sector. In general relativity, this transmission can 
occur only through higher-order processes and is usually suppressed. Though, this trans- 
mission might be substantial in a yet to be found ultra-violet (UV) completion of general 
relativity, or the unified theory may allow the parity violation itself in the gravity sector. 
For example, such phenomena indeed occur in superstring theories. Thus, experimental 
or observational detections of parity violation in the gravity sector certainly provide us 
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with direct information concerning the UV completion of general relativity or the ultimate 
unified theory. 



The most promising route to this end would be to explore parity violation in relics 
of inflation in the early universe. In particular, measurements of parity violation in the 
primordial gravitational waves produced during inflation are expected to bring us informa- 
tion about the Planck scale physics. From this standpoint, it has been intensively discussed 
how to observe the parity violation in the power spectrum directly using the laser interfer- 
ometer pi ^, ||] and indirectly using TB correlation in the cosmic microwave background 
(CMB) [Q, |5|. Remarkably, some speculative gravitational theories with a parity violating 
term produce significant parity violation in the power spectrum of primordial gravitational 
waves |6|, [JJ. Unfortunately, in the conventional inflationary scenario a parity violating 
term leads to only very small amount of circular polarization in the power spectrum of 

primordial gravitational waves i, i, 0, |n], m, m, 0, u, m- 



In principle, we can also seek parity violation in higher-order correlation functions [17, 



18]. Since the power spectrum and the higher-order statistics are often sensitive to differ- 
ent kinds of parity violating interactions, they may have different correlations with other 
revealing statistical features. For example, in a recent paper of Maldacena and Pimentel, 
they discussed graviton non-gaussianity in the exact de Sitter spacetime and found that 
the pattern of parity violating bispectrum is severely constrained by the conformal invari- 
ance jOj. One might then be tempted to deduce that their result implies the unsuppressed 
parity-violation in the bispectrum. We argue that this is not immediately obvious since 
they have not calculated observable quantities explicitly. 

In this paper, we develop a new tool to analyze graviton correlation functions and 
show that parity violation does not show up in the case of the exact de Sitter spacetime. 
We will make clear the origin of the apparent discrepancy between our result and that 



by Maldacena and Pimentel IS]. Furthermore, by extending the analysis to the slow-roll 
inflation, we find that parity violation in the bispectrum does not vanish and its magnitude 
is proportional to the slow roll parameter. Although the parity violation we found is small 
in the case of slow-roll inflationary scenario, it can become significant in the non-standard 
inflations such as Dirac-Born-Infeld inflation [2C, 21]. Much larger parity violation is also 



expected to occur in the recently developed effective field theory approaches 22, 23, 24] to 
which our formulation can be also applicable. 

The organization of the present paper is as follows. In section ||, we present a useful 
formalism for evaluating non-gaussianity in the helicity basis. In section ^ we calculate 
graviton bispectrum in the exact de Sitter universe and explicitly show that no parity 
violation arises in the non-gaussianity. In section [2], we explain how to calculate graviton 
bispectrum in the case of slow-roll inflation. In particular, we demonstrate that parity 
violation in the non-gaussianity exists and is proportional to the slow roll parameter. We 
also discuss observability of parity violation in the conventional inflation and more general 
theories. The final section is devoted to conclusion. In Appendix |A|, we give details of our 
formalism. In Appendix |B], we present useful formulas for polarization tensors. 
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2. A Formalism for Graviton Non-gaussianity 



In this section, with the aid of the helicity basis, we present a useful method to evaluate 
graviton non-gaussianity generated by a parity-violating Weyl cubic term. 

Let us start with the Friedmann-Lemaitle-Robertson- Walker (FLRW) metric 



ds 2 = a 2 (r]) [-di] 2 + 5 i jdx l dx j ] 



(2.1) 



where i,j are indices of the spatial coordinates. No distinction is made between their 
upper and lower indices hereafter for the three-dimensional tensorial quantities. Tensor 
perturbations on this background universe are defined by 



ds 2 = a 2 (rj) [-di] 2 + (Sij + hij) dx l dx j ] , 



(2.2) 



where hij obeys the transverse traceless conditions ha = hijj = 0. The gravitational action 
for the tensor perturbation reads 

Sq = -—7 f dr]d 3 xa 2 
4k 2 J 



h i; j \j hij t khij,k 



(2.3) 



where the prime denotes the differentiation with respect to the conformal time 77 and 
k 2 = 8itG with the Newton constant G. We have two physical degrees of freedom for 
tensor perturbations which can be characterized by the symmetric polarization tensors 
e-^(fc) satisfying 

4f(k) = 0, kjeff(k) = 0, (2.4) 



where k is a comoving wavenumber vector, and s = ± represents the helicity states ±2. 
Namely, they satisfy 



e ijl 



d_ 

dxi 



ikx 



(s), 



ikx 



with k = \k\. It is convenient to adopt the normalization 

e^(k)ef\k) = S ss/ , 



(2.5) 



(2.6) 



where * represents a complex conjugate. Remark that if we choose the phase of each 
polarization tensor appropriately, the following relations hold: 



-ij 



-ij 



eh '(k) 



(2.7) 



These relations will play a crucial role in proving that no parity violation occurs in non- 
gaussianity for the exact de Sitter universe. 

Now, it is straightforward to quantize tensor perturbations. The mode expansion is 
written as 



hij(x,r)) = 2 K J (27r) 3^ E [^(k)u k (r,)a s (k) + e<'\-k)u k {rifa\{-k) 



(2- 
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where the creation and annihilation operators are normalized as 

'a s (k),at(k')] =5 ss ,5(k-k') . (2.9) 
The mode function ut satisfies the evolution equation 

u'l + 2—4 + k 2 u k = . (2.10) 
a 

The canonical commutation relation leads to the normalization condition 

* d d „ 2ik ,.„, 

Uk d^ Uk - Uk d^ Uk = -^- (2 - 11} 

Once a set of mode functions satisfying this normalization is specified, the corresponding 
Fock vacuum is determined by a s (fc)|0) = 0. Then, we can calculate graviton correlation 
functions for each interaction with the help of the standard perturbation technique. 

The purpose of this work is to explore the possibility of parity violation. It is well known 
that the information of gravitational waves is completely encoded in the Weyl tensor W^ v \ p . 
When we restrict ourselves to pure gravity, possible parity-violating interaction terms can 
be found easily. Apparently, the linear term in the Weyl tensor 



d^x^e^W^Xp (2.12) 
vanishes, whereas the quadratic term 

J d^x^ge^W^W^p (2.13) 

is a topological term. 1 Thus, the first two terms in WauXp are irrelevant to the parity 
violation. It follows that the leading term comes from the following cubic action 

Spy = ~b j dr,d 3 x^g-e^ x PW^W a pi S W Xfn s , (2.14) 

which is P- and T-odd. Here, b is a constant with dimension [length] 2 (in c = H = 1 units). 
We shall evaluate the graviton non-gaussianity generated by this term and examine if the 
parity violation emerges in the bispectrum. 

We first establish an efficient formulation to calculate graviton correlation functions 
in a general FLRW universe. In this method, the fact that the helicity decomposition is 
related to the decomposition by a pseudo-duality in the Minkowski spacetime plays an 
important role. To see this, let us consider gravitational waves 7(m)ij i n the Minkowski 
spacetime (see Appendix [A]). In terms of the new variable defined by 

7(M)ij = g ( 7 (M)ij =F ie jkn(M)ik,i) , (2-15) 
which is transverse and symmetric, we can express the cube of the chiral combinations (A.l) 

r (M) 



of the Weyl tensor as (M^ 111 ) 3 = 64(7^, J 3 . This remarkable relation makes computations 



quite simple. 



lr This is related to a Pontryagin number due to R^xp * R^ Xp = W^uXp * W pvXp M 
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When we consider gravitational waves in an expanding universe, however, this 
simple result does not hold anymore. Nevertheless, as is shown in the Appendix A, if we 
define 7^ in terms of the tensor perturbation hij in an expanding universe ( |2.2[ ) as 



(2.16) 



we obtain a useful relation 



1, 



W\ Xp (h) = -W» uXp ( 7 ) 



Minkowski 



(2.17) 



where the Weyl tensor on the right-hand side is understood to take the same functional 
form as that for tensor perturbations in the Minkowski spacetime. Since the algebraic 
structure of the Weyl tensor is the same as that in the Minkowski background, if we define 
by 



1 



lij =F iejkHik,i 



(2.18) 



the relation (VF 1 * 1 ) 3 oc (7 ± ) 3 continues to hold in the FLRW universe. Therefore, we can 
write the parity-violating cubic action (2.14) in a general FLRW background in terms of 
these new perturbation variables 7^ as 



Spy = 8ib / drjd xa 



j3 -5 



(7S)'(7S)'(7£)' - (7£)'(7£)'(7«)' 



(2.19) 



From this, we can read off the interaction Hamiltonian Hpy{rf) as 
8* 



-5 



(fkid i k 2 d i k^5(ki + k 2 + k 3 ) 



7j(^ fc i)) (ijk^faj) (TfciC^fea))'- (l jk {vM) 



(2.20) 



Hj 7 7jj does not represent gravitational waves in 



We should emphasize that 7^ 
Minkowski spacetime but an auxiliary field. From the definition (2.16), we can give an 



explicit mode expansion expression for 7- as: 



d?k 



{2irfl 2 V2k 



Yl [ e if( k ) ( 9 V 7 isk) (au' k ( V )) a s (k) 



+ ef\-k) (d v =f tsk) {au' k {r,)Y a\{-k) 



ikx 



C2-21) 



Practically, this expression can be regarded as the definition of 7^ . Note that it is generally 
impossible to give a closed expression for 7- -. Nevertheless, it does not cause any problem 

because only its time derivative (lifj appears in the interaction Hamiltonian ( 2.2C| ). 

Remarkably, in the exact de Sitter case, we see that 7y has a formal similarity to that in 
Minkowski spacetime and the decomposition 7- corresponds to the helicity decomposition. 
Note that, in general cases, 7^ and 7^ include both helicity components. The main utility 
of our formalism is that a concise expression ( 2. 19| ) considerably simplifies calculations of 
bispectrum in a general expanding universe. 
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3. No Parity Violating Bispectrum in de Sitter Universe 



Employing the formulation developed in the previous section let us compute the bispectrum 
in de Sitter universe. We prove that no parity violation appears in non-gaussianity due to 
the symmetry of de Sitter universe. 

The mode function in a de Sitter background reads 

H 



u k = -(l + ik V ) e~ ikr > 



k 

where H is a constant Hubble parameter. Letting 

d 3 k 



(3.1) 



(27T) 



ikx 



(3.2) 



and substituting the above mode function into (|2.2l|) , we obtain the following formula 



,(+) 



*(-), 



•J 



and 



2i^(2vr) 3 / 2 ^ r 



(3.3) 



(3.4) 



Here, the operator 7^ annihilates a helicity +2 graviton and creates a helicity —2 graviton. 
While, the operator 7^ does the opposite. In the de Sitter background it is easy to integrate 
these equations. In the asymptotic limit r] = 0, we have 



tJ(o,p) 

7,j(0,p) 



-2 K (2^) 3 /\/^ 
y ; V 2 



-2k(2tt) 



3/2 P 



2 L 



e W(p)a + (p) + e*H(-p)al(-p) 
e H(p)a_(p) + e *W(-p)at(-p)' 



Thus, we arrive at the simple relation 

hij(0,p) 



H 



7+(0,p) + 7 n(0,p) 



It is worth noting here the following property 



(3.5a) 
(3.5b) 

(3.6) 
(3.7) 



This traces back to the property of helicity basis (|2.7|) and plays a key role in proving no 
parity violation in the exactly de Sitter universe. 

In the calculations of bispectrum based on the in-in formalism, we need 



(OhtMp) (75fo,*0)' |0> = 2m 2 (27r)W„ (p)<5(fc + p) e *» 



(0| GrSfa.*)) 7j(0,p)|0) = -2^ 2 (2vr) 3 A; 2 n^(p)<)-(fc +P)e 
{0| 7 g(0,p)(7^(^fc))'|0) = (0|( 7fe ± i (r ? ,fc)) / 75(0,p)|0) = 0, 



-ikr) 



(3.8) 
(3.9) 
(3.10) 



Since Tiij and 7^ are related via the differential equation ( |2.1q ), 7^ at 77 = can take arbitrary value. 
Hence, Eq. (3.5) should be regarded simply as the definition of 7^(0, p). 
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where we have introduced the projection operators 



n 



± 

ij,kl 



(p) = e^\p)el ( r'(p) . 



*(±), 



(3.11) 



Given these ingredients, we can compute the bispectrum as 
(74(0,P 1 )74(0,P 2 )74(0,P 3 )) 

d V (0\ \HMV), 74(0,Pi)74(0,P 2 )74(0,P 3 )1 |0> 



\3JJJj 



5! 



(pi +P2+ P3) 6 



where we have used the formula 

i-O 



5! 



{pi +P2+P3) 



(3.13) 



The other non-zero contribution is given by 
(7^ (0, Pi)% j2 (0, P 2 )%j 3 (0, P 3 )) 

= i [ d V (0\ [F PV (r/), 7^ 1 (0,Pi)7~ i2 (0,P 2 )7^ 3 ( 'P3)l |0> 
= -3Mibn*H 5 (2vr) 3 P ^(Pi + p 2 + p 3 ) fa + 1 + ^ 

X n i Ul ,fci(Pl) n i2j2,im^ 2 ^ n j3i3,mA:^3) + ^^(Pl)^ ja.ZmV/^^jg.mfcV 

It can be also verified that the mixed parts (7 + 7 + 7~) and (7 + 7 + 7 _ ) vanish, as they should 
due to the symmetry of de Sitter spacetime. 



M (3-14) 



These results are in agreement with those obtained by Maldacena and Pimentel [19 



They argued that this is determined by the conformal symmetry of de Sitter spacetime. It is 
worthwhile to emphasize that neither ((7 + ) 3 ) nor ((7~) 3 ) themselves are direct observables. 

Now, we need to clarify if the parity violation is observable or not in the non-gaussianity. 
Let us define the right-handed and left-handed circular polarizations by 



h 



R 



hi a 6 



*(+) 



, h L :=h tj $ \ (3.15) 

respectively. A possible observable quantity in which a parity violation is encoded is their 
difference ((h R ) 3 } — ((h L ) 3 }. Here we have 

(h R > L (0, Pl )h R > L (0,p 2 )h R > L (0,p 3 )) 
H 3 

— P . v . ' I TJ-, \i . 

1232 



pIpIpY 113 



e hti (Pi) e S? (PaXfJ (Pa) 



x <(74(0,Pi) +7^ 1 (0,Pi))(7i- 2 (0,P 2 ) +7^- 2 (0,P 2 )) (7^3(0, P 3 ) + 7^ 3 (0,P 3 ))) 



P1P2P3 



x ( (% h {^Pi)%j 2 (O.PajT^aCO, Pa)) + (7^(0^1)7^(0^2)7^3(0^3)) ) • (3-16) 
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For explicit formula, it is convenient to define a function (see Appendix ^ for derivation of 
the last equality) 

F( P1 ,P2, Ps) := 4S + VKi +) (P 2 )e£W 



,*(-), 



{Pl +P2+ P3f {Pl +P2- P3) {P2 + P3 ~ Pl) {P3 +Pl- P2) 

^P\P\P\ 

Then, we obtain the bispectrum of right handed circular polarized modes 
(^(0,^)^(0,^2)^(0,^)) 



*(-)/ 



.(3.17) 



384ib(2Tr) 3 K 6 H 8 F( Pl , p 2 , p 3 ) 



5! 



5! 



(Pl +P2+ P3) 6 (Pl +P2+ P3) 



0. 



HP1+P2+P3) 

(3.18) 



and the bispectrum of left-handed circular polarized modes 

(h L (0, Pl )h L (0,p 2 )h L (0,p 3 )) 
= 384ib(2ir) 3 K 6 H 8 F( Pl , p 2 , p 3 ) 
= . 



5! 



5! 



(P1+P2+P3) 6 (P1+P2+P3) 6 . 



8(Pi +P2 + P3) 
(3.19) 



In the above calculations, the contributions from ( |3.12 ) and ( |3.14| ) have been cancelled out. 
It turns out that there exists no parity violation in a pure de Sitter universe. Since Malda- 
cena and Pimentel proved that *WW 2 is the only way to break parity in the bispectrum, 
we have established that no parity violating bispectrum exists in exact de Sitter spacetime. 



4. Parity Violation during Slow Roll Inflation 

In this section, we extend calculations in the previous section to a slow-rolling inflationary 
universe. The strategy is as follows. We already know there is no parity violation in the 
pure de Sitter case. Hence, when we expand every term up to the first order in the slow 
roll parameter, the sources of parity violation stem from the following three parts: 



(i) change in the asymptotic mode function 

(ii) change of 7^ 

(iii) change of the cosmic expansion 



In what follows we shall discuss these contributions separately. We finally put each con- 
tribution together to conclude that parity violation appears in the first-order in slow roll. 
We will see that the contribution (i) is higher-order in slow-roll, i.e., it is insignificant. 
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4.1 Change in the asymptotic mode function 

In the slow-roll stage, the universe undergoes the following evolution 

ds 2 = -dt 2 + a 2 {t)5 ij dx i dx j = a 2 (ri) [-drf + <%cteW] . (4.1) 

Though the explicit functional form of the scale factor is sensitive to the inflaton potential 
and /or kinetic term, the scale factor takes a simple form in the leading order of the slow-roll 
parameter. To see this, let us define the slow roll parameter e as 3 



H H' __ a a 



" ~ m ~ aH* ' H ~ a~ a 2 ' (42) 

where the dot denotes the derivative with respect to the cosmic time t. In the leading order 
of the slow-roll approximation this slow-roll parameter e can be regarded as a constant. 
Under this condition by integrating 

(■s)'-'- 1 - (43) 

which is equivalent to the definition of e, we obtain 

a(rj) = {-H^y 1 '^ = {-H^fl 2 - V , (4.4) 

where H* is a constant of integration, and we have defined 

3 e 3 , A 

z^ = -H ~- + e. (4.5) 

2 1-e 2 V ; 



The equation of motion for gravitational waves in this background can be written as 

2 

(1 -e)r, 



By solving this, we find that the mode function which satisfies the normalization condi- 
tion ( 2.11] ) and approaches the Bunch-Davis type mode in the rj — > -co limit is given 
by 

where is the Hankel function of the first kind. In the asymptotic limit —rj — > 0, this 
function freezes out and approaches to the constant value 



u k (v) = \e^ 2 



l + e{ 2- 7 -log2 + log-^ 



(4.8) 



where 7 is the Euler constant. Here, we used 



H v ] {z) ~ -- r (^) (7;) V y as |z| < 1 for Re[u] > . (4.9) 

7T V 2 / 



3 This is slightly abuse of nomenclature since the inflaton may not be slowly rolling the potential. The 
parameter e simply measures the departure from constancy of H . 
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and T'(z) = ip(z)T(z), where ip(z) is poly-Gamma function. If we absorb the phase into 
the definition of the ladder operators as a^(k) — > e~ me / 2 a^{k), we find that a relation 
between hij and 7^ at 77 = is the same as that in de Sitter case except for a real-valued 
numerical factor: 

" " ' - (4.10) 



where 



C = 1 + e (2 - 7 - log2 + log 



(4.11) 



As it turns out, this term does not contribute to the parity- violation in the first-order in 
slow-roll parameter. 

4.2 Change of 7^ 

Next, let us consider the change in the expression for the interaction Hamiltonian in terms 
of the creation-annihilation operators. In the lowest order in the slow roll parameter e, 
this change is caused by the change of 7^ from that in the pure de Sitter case in the linear 
order in e. In this section, we denote 7^ in the slow roll inflation background by 7^ to 
distinguish it from the corresponding quantity in the pure de Sitter background, which is 
denoted by 7^ without the bar. 
Now we have 

[d v - isk] (au' k (r,)) =fcVW2-W4y| Q - v\ (-kr,)' 1 / 2 (-k V ) 

+ (-kT]) 1 / 2 Hl% (-krf) + is (-k V ) 1 ^ 2 {-hi) 



(4.12) 



Inserting v ~ 3/2 + e, we can approximate it as 

[d v -isk] (au k (ri)) 
= ik 2 e ine / 2 



(1 + s)e- ikr > + e \ -!— e - ik ^ - — (1 + s)e- ikr > - (1 - s)e lkr 'E 1 (2ikri) 
—kr\ 2 



(4.13) 



where 



and we have used the formulas 



dH ( J\z) 



3v 
dH$\z) 



POO -t 

E 1 {z) = J —dt , 



u=l/2 



dv 




v=-l/2 



2 

TTZ L 
2 

TTZ L 




-iE l {-2iz)e~ iz - ^e iz 



-E 1 (-2iz)e- iz -i^e iz 



(4.14) 

(4.15a) 
(4.15b) 
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Absorbing the phase factor e tnt ^ 2 by rescaling, we have 

d 3 k 



7, 



(27r) 3 /2^2fc L 



4 +) (fc) [d v -ik] {au' k ( V ))a + (k) 
+e£>(-k) [d v -ik] (au'MY aU-k) 



+e\ j \k) [d v + ik] {au'M) a-W 
+e<7 ) *(-*0 [d v + ik\ [au' k ( V )Y al(-k) 



ikx 



(4.16) 



7ij ) e ^ii ' 



where (^7ijJ is the expression in the pure de Sitter case given in ( plf ) and (p^). In terms 
of the functions p± defined by 



p + (k,7]) = i—e~ ik ' n - iTre- lkr ' 



-kr\ 



p_(k,rj) = i- 



-kr\ 



- ikT > - 2E 1 (2ikr ] )e ikr > , 



(4.17a) 
(4.17b) 



the correction term Xij can be expressed as 



d 3 k 



(2ir) 3 / 2 V2k 



^\k)p±{k,r l )a + (k)-e^>' f (-k)p T (k,r l ya\(-k) 



(+)*/ 



+4 \k)p T (k, V )a^k) - e£ j *(-fc)/>±(A;, 7?)*aL(-fc) 



(-)*, 



'ikx 



(4.18) 



The tensor xfj describing the deviation of jij away from that in the pure de Sitter case 
appears in the interaction Hamiltonian as 



-ibe I d 3 xa~ 5 24 



lij) (tJ) Xfci-h.c. 

ibe 24 / d 3 kid 3 k2d 3 k 3 5(ki + k 2 + £3) 
(2vr) D J 

7j(^fel)) (tJ(»7>*2)) XuivM) 



lij(v,ki)^j (j jk (v,k 2 )^ X ki {v,k 3 ) 



(4.19) 



where h.c. denotes the hermite conjugation of the preceding term. 

Now, we can calculate the bispectrum in the slow-roll inflation due to the change of 
7^. From the structure of xfj, we see that graviton non-gaussianity consists of the four 
parts: (7 + 7 + 7 + ) , (7~7~7~), (7 + 7 + 7~)> and (7~7~7 + ). We calculate these parts in order 
using the Wightman functions, 

(0\^Mp)xUv,k)\0) = zK 2 (2n) 3 k 2 Uf hkl (p)6(k+p)p* + ( P ,r ] ) , 



(4.20a) 

(0|xS(^> fe )7^(0,p)|0) = -iK 2 (2irj 3 k 2 nl kl (p)8(k + p) P+ (p, V ) , (4.20b) 
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and 



3 l2tt± 



ij,kl y 



-iK 2 (27r) 3 k 2 Ul kl (p)5(k+p)p^p, V ) . 



(0\xUv,kh±(0,p)\0) 

4.2.1 (+ + +) part 

This part is evaluated as 

(7^(0,^)74(0^2)74(0^3)) 

d v (o\ \h pw (v), ri i Mp 1 hZd 'P2hZiMP3)\ lo) 



(4.21a) 
(4.21b) 



mbe^H^fpipip^ +p 2 + p 3 ) 



+ 5 permutations , 
where A + is the integral 

A+ = 



(4.22) 



drjrj e 



5 JOi+P2>7 „* 



P*+(P3,V) 



-00 




di]7] e 

1 

4! 



—%- 



+ me 



5! 



iir- 



P3(P1 +P2 +P3) 5 (Pi +P2 +P3) 6 ' 



(4.23) 



In contrast to the exact de Sitter case, this integral gives rise to an imaginary part. This 
produces a crucial difference in the final result regarding the parity violation. Here, we 
should stress that both of the circular polarization modes contribute to the (+ + +)-part 
of bispectrum. 

4.2.2 ( ) part 

This part can be obtained from (7 + 7 + 7 + ) by complex conjugation and the flip of the 
directions of each momentum p^. 

(°> (°> (0, P 3 )> 

r ° r 1 

d V (0\ H Py (v), 74(0,Pi)74(0,P 2 )74(0,P 3 ) lo) 



= -mbe^Ht{2TTYpip z 2 pi5{ Pl +p 2 +p 3 ) 



+ 5 permutations . 



(4.24) 



Here, again, we see that both polarization modes make nonvanishing contributions to the 
bispectrum. 
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4.2.3 ( h) part 

In contrast to the pure de Sitter case, we have the cross contributions 



{% lh (0, Pi)% 2h (0, P 2 )%j 3 (°> Pa)) 
i 



f d V (0\ \Hpv{v), 7^(0,^)7^(0,^2)74(0,^3)1 10) 

J — oo 

-mb€K & Hl{2Kfp{plpl5{ Pl + P2 + pa) 

; ^ l! w(Pl) n l " 2J2 ,/m( P 2)n+ j3)mjfc (p 3 )^- + n J 4,W^l)^2,/m( P 2)n i " 3J 3 imfe (P3)^- 

+ {(nii) fe/ 2 )} • (4.25) 



Here, A_ is the integral 



d W 5 e i(pi+P2) Vl(P3,r/) 



— oo 




d-qrfe 

i 

4! 



5 p *(pi+P2)»7 



-i e ip:iV - 2E 1 (-2ip 3 r ] )e 



rp:vn 



-P3V 



P3(Pl +P2 +P3) 



2K, 



where 



I 5„i(pi+P2-P3)»? 



dr]r]°e 



/ — (it . 



(4.26) 



(4.27) 



In this expression, let us set pi + p 2 — p 3 = E and regard £ as a constant independent of 
P3. Differentiation of the corresponding expression gives 



_d_ 

dps 



K 5 



- _1 f° 

PS J-o< 



drjrj e 



5 j(2p 3 +E)r? 



1 5! 

~p 3 {E + 2p 3 f 

d 5 1 1 
dE E p~ 3 ~(E + 2p 3 ) 

d 5 d 1 



log 



dE 5 dp 3 E *\E + 2 P3 



2p 3 



(4.28) 



Since the integral must vanish in the limit P3 — > too, this equation can be integrated to 
yield 

d 5 ( 1 , ( 2p 3 



"F lo g 



dE 5 \E~° \E + 2p 3/ 
Thus, it turned out that A_ contains no imaginary part in contrast to 



(4.29) 
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4.2.4 (+ + -) part 

In parallel with the previous one, this part of the bispectrum reads 



<7,tji (°> PihZi2 (°> P^km (°> P3)> 

= i / d77<0| [i?P V (r/), 7i- 1 (0,Pi)7i 2 (0,P2)7^,(0,P 3 )l |0) 

J — oo 

= 96ibe K 6 H!(2ir) 3 pjp 2 2 p 2 3 6( Pl + p 2 + p 3 ) 

+ {(«ljl) ^ («2J2)} • 

4.3 Change of the cosmic expansion 

The final contribution comes from the imaginary part of the integral 



+ , 

izjz,mk 



P 3 M- 
(4.30) 



/ (ir ? (-r ? ) 5+5e e- i ( pi+P2+Ps )'' 

J — oo 



T(6 + 5e) 



6+5e 



H) 6+5e (pi + P2 + P 3 ) 
5! / 57T . \ = 

~(P1+P2+P3) 6 V + T 7 = "(P1+P2+P3) 6 



5! 



-eB, 
(4.31) 



where the real part of 0{e) has been neglected in the second line since it has nothing to do 
with the parity violation. In terms of this integral, (7 + 7 + 7 + ) and (7~7~7~) are expressed 
as 

<7i- 1 (0,Pi)7i- a (0,P 2 )74(0,P 3 )) 

= ZMibe^Hl{2Kfplplpl5( Pl + p 2 + p 3 ) 



U idiM^ U tj2M^ U tj 3 ,mk(P3) B * + n i 1 j 1 ,fcz(Pl) n i 2 i 2 ,Z m (P2)n i3i3imfc (p 3 ) 



B 



(4.32) 



and 



<7 ilJ - 1 (0, Pi)% 2j2 (0, P 2 )7i 3j 3 (0, p 3 )) 
= -384ibeK 6 HU27T) 3 p 2 1 p 2 2 p 2 3 8( Pl +p 2 +p 3 ) 



U h ju kl(Pl) U i 2j2 M( P ^ U i333,mM 

(4.33) 

From these expressions, it is found that the imaginary part of B fails to cancel the contri- 
butions to the non-gaussianity from (7 + 7 + 7 + ) and (7~7~7~)- 

4.4 Parity Violation 

Now, we are in a position to discuss parity violation in a slow-roll inflationary universe by 
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assembling results obtained above. To see the parity violation, we are required to evaluate 
(h R ' L (0, Pl )h R > L (Q,p 2 )h R > L (0,p 3 )) 

P1P2P3 



x ((7 < t ix (0,Pi)7i- a (0,P2)7j i3 (0,P3)> + (74(0,P 1 )7^ 2 (0,P 2 )7^ 3 ( 'P3))) 



(4.34) 



In the previous subsections, we found that the three-point correlators of 7^ are of order 
e. It follows that the change from the overall constant factor C becomes higher-order in e. 
Thus we set C = 1 henceforth. From Eqs. fl422p , (g3|), flOp , and (|4~33|) , we obtain the 
final result 

(^(0,p 1 )^(0,p 2 )^(0,p 3 )) 

= -32i6e(2vr) 3 K 6 J ff, 8 F(p 1 , p 2 , p 3 ) [34+ - 3^; - 2B + 2B*} 6( Pl +p 2 + p 3 ) 
+5 permutations 

6! 

= -64(27r) 4 e& K 6 F^( Pl + p 2 +P 3 )F(p 1 , p 2 , p 3 ) , - '— , G , (4.35) 

{Pi +P2+ P3) 

and 

(h L (0, Pl )h L (0,p 2 )h L (0,p 3 )) 

= 32ibe(2Tr) 3 K 6 HfF( Pl , p 2 , p 3 ) [3A+ - 3A; - 25 + 25*] <5( Pl +p 2 + p 3 ) 
+5 permutations 

_6!_ 

(pi +P2 +P3) 

where F(pi,p 2 ,p 3 ) has been defined in Eq. ( 3.17] ). Apparently, parity violation shows 
up in the bispectrum and its magnitude is proportional to the slow roll parameter. It is 
interesting to observe that the bispectrum of curvature perturbations is also proportional 
to the slow-roll parameter in the conventional single inflationary scenario [17]. 

It should be noted that no parity violation occurs in the mixed parts. Introducing a 
function 

J(pi,P2, P3 ) := et\k 1 )et\k2)e^\k 3 ) = e^fex)^ (fc 2 )e^ (fe 3 ) , (4.37) 

a simple algebra shows that 

(h L (0, Pl )h L (0, P2 )h R (0, P3 )) 

= m2(2^febK & HlJ( Pl , P2 , P3 ) [A_ - A*_] 5( Pl + p 2 + p 3 ) . (4.38) 

Since A_ is real, there exists no parity violating contribution. The same conclusion applies 
to the bispectrum (h R h R h L ). 

It should be stressed that the parity violation can be observed in the CMB. In fact, 
three-point correlators (TTB) , {TEB) , (EEB) become non-zero in contrast to the parity 
conserving cases. In the conventional slow-roll inflationary scenario, the amplitude might 
be too small to be detected in near future. However, in the non-conventional scenarios, we 
might have much larger parity violation. 



M(2n) 4 ebK 6 H s J( Pl + p 2 +p 3 )F(p 1 , p 2 , p 3 ) — ; - ; , (4.36) 
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5. Conclusion 



We have developed a useful formalism to evaluate graviton correlation functions. As an 
illustrating application of this formalism, we studied parity violation in the early universe 
through non-gaussianity of gravitons. First of all, by calculating the bispectrum we have 
shown that no parity violation arises in the exact de Sitter background. This statement 
may appear to be inconsistent with that by Maldacena and Pimentel. However, the explicit 
results of calculations they gave in their paper are in agreement with ours. The difference 
came from the interpretation. They have calculated ((7 ± ) 3 ) and concluded there is a 
parity violation. This is not the case since the actual observables are ((h R ) 3 } — ((h L ) 3 ), 
etc. These vanish in the purely de Sitter universe. The situation is different when the 
spacetime departs from the exact de Sitter. In slow-roll inflationary case we have found 
parity violation in the graviton bispectrum proportional to the slow roll parameter. 

We also discussed that parity violation in the bispectrum can be observed e.g., in 
the (TTB) correlation in the CMB. It might be also possible to detect the signature of 
parity violation through direct observations of primordial gravitational waves using a space 
interferometer observatory [26]. 

Although we have concentrated on the parity violating Weyl cubic term, the present 
formalism enables us to calculate higher-order parity violating correlation functions rather 
straightforwardly. In fact, it does not require a sophisticated expertise to evaluate the 
parity violating non-gaussianity of the type *WW n . We can also extend our analysis 
to more general inflationary scenario using the effective field theory approach. In those 
cases, we expect large parity violating non-gaussianity. Moreover, we can evaluate parity 
violating contribution due to a non-minimal coupling to the inflation <j> such as a term 
f{4>) * WW [0]. Our formalism would be useful in analyzing these general cases. Work 
along these directions is in progress. 

It is known that exact de Sitter correlation function is related to correlation function 
of stress tensor in the conformal field theory through analytic continuation [17]. Hence, 
the graviton non-gaussianity has an interesting application. Similarly, our results would 
be useful for calculating three-point functions of stress tensor in the non-conformal field 
theory using gravity/field theory correspondence. 
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A. Pseudo Self-duality, Helicity, Parity Violating Action 

The pseudo self-dual and anti-self-dual Weyl tensor are defined by 

W± Xa := W puXa ± i * W^xa, *W^ Xp = l -e^W apXp , (A.l) 

which satisfy 

*W%Xp = TiW± Xp . (A.2) 

The factor "i" arises due to the Lorentzian signature. Since the relation W pu a PW~p lS = 
holds, we get 

e^W^W a ^W XpiS = 1 (w+ a ?W+/W+^ - W- u ^W-/W~r) ■ (A.3) 

Note that we also have a relation 

W^W^W^ = i + W- "^"/W^) • (A.4) 

Let us consider tensor perturbations in the Minkowski spacetime 

ds 2 = -drf + (5ij + 7(M)y) dx l dx j , (A. 5) 

where J(M)ij is a transverse traceless symmetric tensor describing a gravitational wave. 
The on-shell equations of motion are 

-Vi4)tf + V 2 7(M)tf = 0, (A.6) 

where the prime denotes the differentiation with respect to r\. It is convenient to define a 
new variable by 7(M)ij := e jki7(M)ik,i- By multiplying e m jj to both side, we can prove this 
new tensor is symmetric. Then, it is easy to see the transversality of 7(M)ij- Thus 7(M)ij 
is also the transverse traceless symmetric tensor. Using this new tensor, we can define 
helicity eigenstate by 

T(M)ij : = 2 ( 7 (M)y ^ *7(M)tf) ' ( A - 7 ) 

One can then verify the relation 

g ± 

A direct calculation yields 

W°m = ^7(M)i fc (A.9) 

W°jfcZ = -7(M)j[fc,/] = -^Iml^jrn (A. 10) 

W^kim = -<%7(M)fc] m + <W7(M)fc]Z » (A- 11 ) 

where we used on-shell equations for gravitational waves ( |A.6| ) in Minkowski spacetime. It 
is easy to derive the following relations 

4 In appearance the expression of W 3 derived here disagrees with Eq. (2.14) in [[L9|. However, they can 
be shown to be equivalent by using three-dimensional identity 2(7( M )i[jfc,j]7(MMk,i] + 7(M)fe[i,i]7(M)*[i,j]) = 

c / / 

°ji7(M)fc(,m7(M)fc[i,m] • 
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and 



W±^W^ s W±r = 64(7 ( ± M)iJ ) / (7 ( ± M) jfc )'(7 ( ± M) fe J 



W^W^W^ = 27(M)ij [7(M)jkl(M)ki - 3 7(M)ifc7(M)fci 

e^ x PW^ al3 W a p 1& Wx P1 8 = 4 7(M)ij 7(M)jfc7(M)fc i ~ ^l{u)jk^{u)ki 



(A.12) 

(A.13) 
(A.14) 



These equations manifest the symmetry between W 3 and *WW 2 in a flat space. 

Now, we would like to promote the above beautiful relation to FLRW spacetime (|2.2|), 
The Weyl tensor is obtained as 



1 



iz ( ah ik)' 



2n 



W° 



jki 



' h 'm\ 



w j k im = - 

a 



- 6 l[j( ah 'k]m)' + 5 m[j( ah, k]l)' 



(A.15) 
(A.16) 
(A.17) 



where we have used on-shell equations for gravitational waves in FLRW spacetime. Com- 
paring this expression with the one in Minkowski spacetime, we notice that if we make the 
identification 

iij = <j , (A. 18) 

a simple relation 



nvW^n^l^ (A.i9) 

holds irrespective of the expansion history of the universe. Therefore, with a definition 

itj = 2 Hj =F iejkHik,i) 

we obtain 

(A.21) 



(A.20) 



W^(h)W^(h)W^(h) = ^yjfr? 
Thus, we can deduce the action in the main text. 
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B. Polarization tensors 



We fix the representation of polarization tensors ^7j. Because of momentum conservation 
k\ + &2 + &3 = 0, we can make all hi lying on the (x, ?/)-plane without losing any generality. 
It follows that a triangle can be constructed as 



k>\ = k\ (1,0,0) , &2 = &2 (cos 9, sin#, 0) , k% = k% (cos 4>, sin <f>, 0) 



(B.l) 



where 



cos I 



h 2 - k 2 - k 2 
^3 fc l ^2 



sm ( 



2kxk 2 



cos < 



'OJ.2J.2 I 01.2J.2 i OJ.2J.2 _ L.4 _ j.4 _ u4 

^jrv-^rx/Q \ ■""'2 3 3 1 1 2 3 



k 2 - k 2 - k 2 

/i 2 ^3 



A 



2k 3 k x 



sm ( 



(B.2) 
(B.3) 
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Note that < 9 < it and n < (ft < 2tt. Using the above representation, we can fix the 
polarization basis. For k\, we have a simple expression 



,(«0 



(kl) 




(B.4) 



where s± distinguishes the helicity. To obtain the polarization for k 2 , we need to rotate 
the above one by 9 as 



e^>(k 2 ) = 

Similarly, we obtain the polarization tensor for fc 3 as 



sin (7 
- sin 9 cos I 
y —is 2 sin9 



sin v cos ( 

2/3 



cos 
IS 2 cos 




(B.5) 



>3) 



(**) 



/ 


sin 2 <j) 


— sin (j) cos </) 


— is 3 sin 




— sin 4> cos 4> 


COS 2 </> 


is 3 cos <; 


V 


—is 3 sin <f) 


ZS 3 COS 


-1 



(B.6) 



From these expressions we can calculate F(k\, k 2 , fe 3 ) and J(k\, k 2 , fe 3 ) defined respectively 
by Q3.17 ) and ( 4.37 ). The resultant expressions are given by 



F(k u k 2 ,k 3 ) 



and 



J(ki,k 2 ,k 3 ) 



-(1 - cos0)(l - cos0)(l - cos(0 - 6)) 

{(h + k 2 f - kl) {(k 2 + fc 3 ) 2 - k\) {(k 3 + fcQ 2 - kl) 
64fe 2 feffef 

(ki + k 2 + k 3 f (ki + k 2 - feg) (fc 2 + fc 3 ~ fei) (fc 3 + fei - k 2 ) 

64fe 2 fe|fe 2 



- - (1 - cos B) (1 + cos <j>) (1 + cos(<£ - 0)) 
8 

(fei + fc 2 + fc 3 ) (fci + fc 2 - fc 3 ) 3 (fc 2 + fc 3 ~ fei) (fc 3 + fei ~ fc 2 ) 

64fc 2 A|£; 2 



(B.7) 



(B. 
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